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The subject of event-by-event fluctuations of charged particles produced in high-energy heavy-ion collisions has 
recently gained considerable attention [|l], |[ |[ ^, |^. The focus has mainly been in the identification of a measure 
that can differentiate a quark-gluon plasma (QGP) from a hadron gas (HG). No consideration has been given to the 
fluctuations introduced by the phase transition (PT) itself. If the FT corresponds to the second-order chiral transition 
that has the usual 0(4) symmetry of two massless quarks, it has been argued that long-range correlation typical of 
critical phenomena is unlikely to occur because the only scale in the problem is characterized by the realistic pion mass 
that is roughly the same as the critical temperature As a result one does not expect large fluctuations, especially 
in the charge sector. On the other hand, if the chiral transition takes place rapidly far from thermal equilibrium, then 
the disorientation of the isovcctor order parameter can lead to long- wavelength modes and large fluctuations 0, |[ D ■ 
The usual signature for such disoriented chiral condensates (DCC) is the large fluctuation in the neutral-to-charge 
ratio of the produced particles |ll|, In this paper we study the nature of the fluctuations within the 

charge sector and identify measures that can be computed explicitly as well as being amenable to direct experimental 
verification. They do not depend on the existence of long-range correlations that may be suppressed by the rapid 
expansion in a heavy-ion collision. 

The usual starting point of a consideration of DCC is the linear a model, for which the potential is 



4 



,2^2 



Ha, 



(1) 



where the chiral fields are represented by a vector $ = (fi^) in 0(4) space. The parameters A,w, and H are 
determined by the masses m.,:,ma and the pion decay constant f-^. At temperature T far below the critical Tc, the 
normal vacuum is characterized by (cr) ^ 0, and (tt) = 0. At T well above there is approximate 0(4) symmetry 
for which ($) = 0. If the QCD plasma is close to thermal equilibrium as it cools from above Tc to below, the PT 
results in ($) becoming nonzero in the a direction with no large fluctuations expected in the angular deviation from a. 
However, in the quench scenario far from thermal equilibrium [Q, the plasma loses touch with the vacuum orientation 
and (tt) becomes nonzero along arbitrary directions in isospin space in different spatial regions, thereby generating 
large fluctuations in the charges of the pions produced. What we seek are the signatures of those fluctuations that 
are independent of the details of the theory, specifically, the parameters governing the chiral transition. We show 
the existence of a numerical index v that can serve as a signature of DCC. As an alternative to the D measure that 
has been suggested in [0, ^ , we consider another measure B that can clearly distinguish the different types of charge 
fluctuations. 

It is convenient to start with the coherent-state representation for the statistical fluctuations since the multiplicity 
distribution of a pure coherent state \a) is Poissonian, i.e., |(n|Q:)|^ — P^^, with average multiplicity 
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where the property a{z)\a) = a{z)\a) has been used. We generalize this formalism by incorporating isospin and treat 
(j) as the eigenvalues of the isovector annihilation operators d{z) Q. The total average density of hadrons (assumed 
to be pions only) is then 
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The spatial coordinate z can be regarded as (pseudo) rapidity in heavy- ion collisions, as in p^ , but at this point that 
identiflcation is unnecessary. Our theoretical results, expressed below by Eqs. (16), (17) and (21), are independent of 
what z is exactly. For experimental analysis of the data, extensive discussion will be given below. 

The application of coherent states and the Ginzburg-Landau formalism |l5| to multiparticle production was con- 
sidered many years ago 0. We use the same approach here, as in ||l^, |l^, to study the thermal fluctuations in the 
chiral transition to DCC. We take the Ginzburg-Landau free energy to be 
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where only the isosymmetric part of the potential in Eq. (|l|) is adapted here. The derivative term in (4) is neglected 
here, since our earlier studies of isoscalar 4){z) indicate that the inclusion of d(j){z)/dz term in the free energy leads 
to negligible effect on the scaling result jlj, Q. When T is lowered below Tc, a becomes negative, while b remains 



positive, and the system makes a transition to the hadron phase whose density fluctuates around 
hadronic multiplicity distribution is then given by 

P(n_|_, n_, no) = ■Z^ "'^ / 'Dcfi P'^{n^,n^,no, 
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Since these quantities are meaningful only in the hadron phase, we shall in the following consider only the situation 
where a < 0. 

The central aim of our proposed analysis is to find a measure of the charge fluctuations due to a chiral transition 
to arbitrary (t?) directions with the property that the measure is independent of the details, more specifically, the 
parameters a, b, and S in Eq. (Q). We shall assume that 6 can be arranged to be small so that (f can be regarded as 
constant inside 6. We then have 



F[(j)] = S a (j) +b (j) , (9) 

where the value of each of |0j| is allowed to vary throughout the complex plane in the functional integrals in Eqs. (|^) 
and d). 

In our search for a quantity that is independent of S, a, and 6, we first consider the bivariate factorial moments 
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for integer values of qi and q2. If P(ri_|_, n_, np) were the statistical distribution given in Eqs. @) and §), then (0) 
would yield fqi^q2 = /o\/i ^'^^ 1^ ^^"^ 92- Deviation from this trivial result on account of Eqs. ^ and (^) is then 
a measure of the effect of the DCC on the charge fluctuations plus other effects to be discussed below. Making the 
appropriate substitutions we obtain 
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Changing the integration variables to the set s — — \(t)+\^ + and u = so that only u is integrated 

from to 00, we obtain a significantly simplified, closed form 
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where _B(m, n) is the Euler-beta function, and 



X — \a\ y 5 /h 



The function Jp{x), defined by 
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can be related to the parabolic cylinder function, but is straightforwardly computable. 

Although /gi,g2 has complicated dependence on S, a and b, different (91,92) moments have the same dependence 
if gi + 92 = 9 is the same. Furthermore, the factor ((5/6)('?i+'?2)/2 cancelled for the normalized factorial moments, 
which we define as 
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Evidently, Fq-^^q^{x) is a function of x only. Its dependence on x is shown in Fig. 1 in a log-log plot. No scaling 
behavior can be seen. However, if we plot InFq-^^q^ vs ln_F2,2 as in Fig. 2, we find a substantial region in which the 
relationship is linear. In that linear region we can write 

F,,,,, cxF^^'^"", (16) 

which is a behavior that is independent of x. More explicitly, we can determine Pq-^^q^ by straightline fits of the curves 
in Fig. 2 in the region 0.2 < lni^2,2 < 0.9. The result yields /3g^ as a function of only the sum, 91 + 92 = as can 
be seen directly from Eq. (^). That dependence of (3q on q is shown in Fig. 3. Apart from the point for g = 2, it is a 
linear dependence of In (3q on ln{q — 1). Thus we can write 

/3, cx(g-l)^ z/=1.29. (17) 

This index v is independent of x, and therefore of 6, a and b (so long as a < 0). The scaling behaviors ( p^ and ( [T7| ) 
summarize the properties of charge fluctuations in a chiral transition to DCC and culminate in a numerical index ly 
that characterizes the phenomenon. 

To verify the above behavior experimentally, it is necessary to vary x, which is the implicit variable in Eq. (p^). 
Since a and b in Eq. (^3|) are not subject to experimental control, only S can be varied. A more extensive discussion 
of the experimental cuts that are optimal for detecting the signal is postponed until another prediction independent 
of X is presented. 

A few remarks should first be made regarding the sensitivity of the above theoretical result to statistical fluctuations 
in the experimental background. It was shown by Bialas and Peschanski p9[[_ that the factorial moments filter out 
the statistical fluctuations represented by the Poissonian distribution in EqT^). However, if the DCC is produced 
in a background of conventionally produced hadrons, we have to consider an additional contribution to the mean 
multiphcity. Thus we replace on the RHS of (||) by |(/>ip -h/S, where S denotes the mean density of the statistical 
background. If there are statistical fluctuations in the orientations of the fields after the chiral transition, we shall 
let 5* represent them also. We ask whether our result is sensitive to a small perturbation by S. Clearly, if S were 
large, it would be hard to find the signature of DCC in the presence of a large non-DCC hadronization process. If 
S is small, we can carry the above theoretical calculation to first order in 5', and find that a term must be added to 
Pqi,q2 ™ (|l5|), which is proportional to SVbd. At small S the effect of S is expected to be negligible. 

There is, however, a limitation on how small S should be allowed to be. That can be seen from Eq.(^) where 
the exponential damping term becomes ineffective at small 6. Without letting |0p become large enough to require 
higher-order terms in the Ginzburg-Landau free energy, a rough lower bound on 6 is derived in Rcf. |l^ and is equally 
applicable here; it is 5 > x^b/a?, where Xq = V41n2 = 1.67. That bound translates to S^/bS > {Sb/\a\)xQ. Thus 
the sensitivity to the statistical background depends on its strength relative to that of the dynamical strength |a|/6, 
as is reasonable. Our conclusion is therefore that unless the mean multiplicity of the statistical background is small 
compared to that of the DCC production, our approach of searching for scaling behavior would not be successful. 

To put the curves in Fig. 2 to experimental test, one inevitably has to deal with the inaccuracies in the determination 
of and in particular with F2.2. If the error bars are large, the effectiveness of this approach is clearly limited. 

Even if they are not large, the data points may admit a larger scaling region than what we have consider in Fig. 2 in 
the derivation of = 1.29. In fact, our lower bound on 5 mentioned above implies that the curves in Fig.l are reliable 
only for x < xq, or — Inx < —0.51. It translates to a restriction on the range of ^2^2 to the region lni^2,2 < 0.65 in 
Fig. 2. Thus fitting the curves in that region by straight lines with allowance for errors we find that it is necessary 
to modify the value of v to 

1.42 ±0.13. (18) 

It is clear that we can no longer claim strict x- independence in our result, when the various complications discussed 
above are taken into account. However, despite the limitations arising from theoretical considerations, we feel that 
the proposal for the experimental measurement of Fq-^ q^ may nevertheless yield interesting insight into the formation 
of DCC. 

To reduce the sensitivity to the details in the scaling analysis, we now consider a global measure that makes contact 
with the observables proposed in ||l|, |^, ^, |^, ||]. Let us first list the simple identities: 

/i,o = ("+) , /o,i = {n-) , = {n+n_) , 

/2,o = ("+(^+ - 1)) , h2 = {n-{n_-l)). (19) 

In terms of the usual definitions, A'ch = + n_, Q — n+ — n^, and (SX^) = {X"^) — {X)^ , we then have 

(<5Q2) =2(/2,o + /i,o-/i,i), (20) 
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where the symmetry fij — fj^i has been used. Let us now consider the quantity 
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Using Eqs. @ and d^), together with (|T|) we obtain 
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On the other hand, in terms of D 
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It has been argued that D < 4 whether the thermal system is a QGP or a HG so in that scenario we have B < 
whatever (n^ri-) / {Nch} may be. That prediction is distinctively different from 2 in our scenario. Qualitatively 
speaking, the charge fluctuations in chiral transition to DCC are much greater than what can be expected from the 
thermal fluctuations in either the QGP or the HG phases. Thus the measurement of B can decisively determine the 
nature of the FT that the QGP system undergoes. The presence of statistical fluctuations in the background will 
undoubtedly weaken this claim to an extent that is proportional to the strength of such fluctuations, as we have found 
in the study of the scaling behavior. 

In order that our theoretical prediction can be applied to the data, it is important that the data are analyzed in 
the proper way, which is the subject we now address. The factorial moments can reveal scaling behavior if they can 
capture the rare events with large spikes without being overwhelmed by the contributions from average multiplicities 
[p9| . In other words, if we, for simplicity, consider a single factorial moment 



f,^J2 "('^ ~l)---{n-q+l)Pn 



(25) 



we see that fq probes the high end of the distribution P„ with n > q; thus if q ^ (n), only the less frequent events 
with n ^ (n) can contribute to Eq. (|2^). That can be achieved for heavy-ion collisions either by using extremely high 
q, or by making severe cuts in data selection to reduce (n). Neither was done in previous analyses of the nuclear data 
and nothing of interest was found in F2 or F3 in the bulk data. 

To test the power-law behaviors (|l^)-(|l^) and the prediction i? = 2, it is necessary to make cuts in the data to best 
reveal the features of this analysis. To see what cuts to make, we first discuss multiparticle production in heavy-ion 
collisions. If the quark-hadron FT is first-order, then what we have considered here is irrelevant. For a second-order 
chiral transition it is likely that, if the DCCs are to be created at all, they would appear as clusters of hadrons in 
different regions in space and time amidst other patches of hadrons produced as a result of the gradual cooling of 
the QGF that undergoes the normal critical transition. If it is possible to restrict the observation to a very short 
duration in the emission time, then we should see many regions of voids where no particles are produced, separating 
the regions of particles produced either as oriented or disoriented condensates. Unfortunately, a selection in emission 
time is not experimentally feasible. The event multiplicity integrates over all times and smooths out the fluctuations 
in clusters and voids. Thus to detect the DCC by our fluctuation analysis it is necessary to make a severe cut in pT 
and select only the particles emitted into a narrow px window. Such a Apx cut achieves two goals: one, to reduce 
the multiplicities to facilitate the factorial moment analysis discussed above; two, to minimize the overlap of hadron 
emissions at different times. The correlation between pT and the emission time is discussed in Ref. |t2C| , where a 
simulation is done to exhibit the void patterns. 

Now, suppose that in a narrow pT window n+ and ri_ charged particles are collected in in azimuthal angle and 
Ay = 6 in rapidity. What we suggest is that (A^ch) be only about 2 so that Fq-^^q^ can be calculated for A^ch > 9 = 91+92 
up to q about 10, as shown in Fig. 1. Such values of A'ch are large deviations from (Nch), but represent only a tiny 
fraction of the total multiplicity produced in a typical event in heavy-ion collisions. Thus the problem of overlap from 
different emission times is minimized, while the characteristics of charge fluctuations that convey the nature of the 
chiral transition are retained. Of course, (Nch) depends on the sizes of the cuts. But our theoretical analysis based 
on the Ginzburg-Landau formalism indicates that the scaling exponent I3q-^^q2 , the index and the quotient B are 
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all independent of the window sizes, so long as there is enough statistics to render accurate determination of fq^.q2- 
Since the charge fluctuations of DCC are much larger than those among the hadrons produced in the normal critical 
transition, the contribution of the latter to our measure is expected to be negligible. 

Recently, WA98 presented their data that reveal no correlated charge-neutral fluctuations in Pb-Pb collisions 
Unfortunately, they have no information on pT and cannot make the px cut that we regard as essential. The absence 
of any signature of DCC in the bulk events does not preclude the possibility of local condensates of the type we 
consider here. 

It is clear from Eqs. ( pO[ ) and ( ^ ) that B involves order-sum q not greater than 2, whereas the power-law behavior 
of f3q shown in Fig. 3 is for q > 2. Although the former seems simpler to investigate, the smallness of q requires 
(-^ch) ^ 1 in order for the factorial moments to exhibit the fluctuation properties at the far edge of the multiplicity 
distributions. That in turn requires the kinematical cuts to be very severe so that only in rare events can N^h exceed 
q. As a consequence, the statistical errors may be large. For the latter study of f3q, the values of q are larger, so the 
cuts can be less severe to allow (A'ch) to be larger. The experimental errors on Pq then may hopefully be small enough 
to render the determination of i' feasible. Whether q is small or large, our analysis probes large fluctuations from the 
mean. 

Despite the many complications that accompany heavy-ion collisions, our analysis suggests a possible way of flnding 
the DCC. If only the bulk events are examined, it is quite possible that the DCC, even if created, might escape 
detection. The simple consideration done here based on the essentials of chiral transition offers a potentially effective 
tool to initiate a flrst-round experimental exploration. Since the discovery of a DCC would provide a deflnitive 
confirmation of the accepted ideas of chiral dynamics in a dense system, such an exploration is well worth undertaking. 

We wish to acknowledge helpful communications with Drs. Q. H. Zhang and V. Koch. This work was supported, in 
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FIG. 1: -Fqi.ga vs X for various orders of the moments. 
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FIG. 3: Power-law behavior of (3, 



